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Abstract 

Finite element models of the thorax are being developed to assist engineers 
and vehicle safety researchers with the design and validation of counter- 
measures such as advanced restrain systems. Computational models have 
become more refined with increasing geometrical complexity as element size 
decreases. These finite element models can now capture small geometri- 
cal features with an attempt to predict fracture. However, the bone ma- 
terial properties currently available, and in particular the rate sensitivity, 
have been mainly determined from compression tests or tests on long bones. 
There is a need for a new set of material properties for the human rib cor- 
tical bone. With this objective, a new clamping technique was developed 
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to test small bone coupons under tensile loading. This technique allows for 
applying minimal constraints to the coupon during clamping and ensures 
that the main type of loading is tension. Ten coupons were harvested from 
the cortical shell of the sixth and seventh left ribs from three cadavers. The 
coupons were tested to fracture under quasi-static (target strain rate of rate 
of 0.07 %/s) and dynamic loading (target strain rate of 170 %/s). Prior 
to testing, each coupon was imaged with a computed micro-tomograph to 
document the bone microstructure. An optical method was used to deter- 
mine the strain field in the coupon for the quasi-static tests. The rib bone 
coupons were found to be elastic, with brittle fracture. No plastic behavior 
was observed in this test series. The bone coupons were assumed isotropic, 
homogeneous and elastic linear, and the average Young's modulus for the 
quasi-static tests (13.5 GPa) and the failure stresses (quasi-static: 112 MPa, 
dynamic: 124.6 MPa) were in line with published data. Fracture however 
did not always occur in the gage area where the cross-sectional area was the 
smallest, which contradicted the assumption of isotropy and homogeneity. 
The comparison of the results obtained in the current study with published 
results on tibia and femur bone coupons suggests that the effective cross- 
section has an effect on the calculated material properties, and that further 
analysis of the bone microstructure is required to establish the rib cortical 
bone fracture mechanism. 

Keywords: bone, cadaver, rib, strain, porosity, fracture, quasi-static, 
dynamic, micro-computed tomography 



1. Introduction 



Rib fractures are a good indicator of the severity of an impact to the tho- 
rax as the protection to the internal organs such as the lungs and the heart 
is greatly reduced with the increasing number of fractured ribs (abbreviated 
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olds under diverse load conditions. A significant milestone was achieved in 
the characterization of the strength of the thorax by accounting for the geo- 



2 



metrical variations in the rib cage and the rib themselves, and for the effects 



of biological variations such as aging (Berthet et al. 2005 Ito et al. 2009 



Gayzik et al. 


2008 


Kent et al. 


2004) 



geometry that includes variation in the shape of the cross-section along the 



rib axis (Kindig 2009), an increase of the twist from the posterior to the 



anterior aspect ( Mohr et al. , 2007 ) , as well as a non uniform distribution 



of the cortical thickness (Choi and Lee, 2009). In a recent study, Li et al. 



(2010b) investigated the sensitivity of the rib structural response obtained 
from a computational model to (a) the accuracy in the reconstruction of the 
rib (quantified by the mesh density), (b) the cortical thickness distribution 
and (c) the material properties. The response of the rib finite element model 
was found to be little sensitive to the choice of material properties, whereas 
clear trends were observed in the effect of the mesh density and cortical 
thickness distribution. In finite element modeling, the cortical shell is de- 
fined as a continuum made of pure cortical bone ( |Li et al" 2011). However, 
rib cortical bone has voids and is not homogeneous (figure [TJ, and therefore 
the bone material properties documented for homogeneous bone are likely 
to inadequately represent the mechanical behavior of rib cortical bone. 




Figure 1: Cross-section of the left 6 th from the antero-lateral aspect, from micro-computed 
tomography. The arrows indicate voids. The bar scale is 1 mm long. 



Bone tensile material properties have been reported extensively for bovine 
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experimental tests by Kemper et al. (2005) are the only study targeted at 
human rib bones. Delicate "dog bone" coupons (length: 30 mm, width: 
9 mm at the ends, 2.5 mm in the gage area, thickness: 0.5 mm) were ma- 
chined from cortical rib bones, and a slack adaptor system was designed to 
apply a constant strain rate, with a target of 50%/s. Special care was taken 
to prevent the breaking of the samples during handling. However, for more 
than 80 % of the tests (95 out of 117 tests) the bone samples fractured in the 
grip area. For the other tests, the actual fracture location was not reported 
and therefore one can only assume that the fracture occurred in the gage 
area. Dog bone samples have a weak point by design - the gage area where 
the width is much smaller than that at the ends of the sample - so as to 
generate greater tensile stress in the gage area and ensure that fracture is 
caused by tension, at least under the assumption that the sample is homo- 
geneous. The fact that fracture occurred in the grip area so often suggests 
either that the test apparatus caused the fracture, for instance because of 
inadequate clamping force or misalignment of the top and bottom grips, 
or that the hypothesis of homogeneity did not hold. In addition, |Kemper 



et al. (2005) used an extensometer to measure the average strain for the gage 



area that may have interfered with the mechanical response of the samples. 
The questions raised in Kemper et al. (2005) and the paucity of rib bone 
material data highlighted the need for new experimental data. The goal of 
the current study was to perform new tensile tests of rib bone coupons to 
link the fracture characteristics to the coupon bone microstructure. To do 
so, an experimental protocol was developed to machine and test rib cortical 
bone coupons of constant thickness under quasi-static and dynamic load- 
ing. All the samples were imaged prior to testing using a micro-computed 
tomograph. This protocol was evaluated on coupons machined from the left 
sixth and seventh ribs harvested from three post mortem human subjects. 



2. Materials and methods 

2. 1 . Coupon preparation 

The ribs were first cut into five 60-mm long pieces named A (most an- 
terior piece) to E (most posterior piece). Isotropic images similar to figure 
[T] were obtained with a micro-computed tomograph (Scanco Medical Viva 
CT40, Bruttisellen, Switzerland) at a resolution of 30 |a.m/pixel and used to 
determine where the bone coupons could be machined based on the curva- 
ture and the thickness of the cortical shell (Figure[2]-a). Next, the periosteum 
was removed, and the sample was shortened with a small abrasive band 
saw equipped with a diamond blade (Diamond bone band saw, Mar-med 
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inc., Cleveland, OH, USA). The next cuts were performed with a low speed 
diamond saw (IsoMet low speed saw, Buehler Ltd, Lake Bluff, IL, USA, 
equipped with a 15HC Diamond blade) that can execute parallel cuts with 
micrometric accuracy. The sample was glued to a piece of hard foam that 
was then clamped in the sample holder in the desired orientation thanks to 
screws with spear tips (Figure [2j-b) . This step was critical in the machining 
process: the orientation of the bone slab had to be adjusted so that the cut 
performed with the low speed diamond saw would be entirely in the cortical 
shell. A first cut was made to generate the outermost flat surface. From 
there, an other cut was made parallel to the first one to create a flat piece 
of cortical bone (referred to as wafer) that was 0.5 mm thick (Figure [2jc) . 
Next the two holes used to reference the position of the bone coupon (figure 
|3|) were drilled in the wafer. 




Figure 2: Machining of bone coupons. 
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Figure 3: Dimensions of the bone coupons. The 2 holes are used only for references for 
machining and testing, and do not carry any load. 

The next step consisted in machining the wafer to give it its final shape. 
This happened in two successive actions. First, the wafer was shortened 
to the desired length, and second, the narrow section of the coupon was 
machined. These two actions utilized the same principle: the wafer was 
sandwiched between two pieces of plastic (Figure [2]-d), and attached thanks 
to two screws to a template made of brass. A router equipped with a bearing 
flush trim bit was used to trim the bone wafer to the desired shape: the 
bearing was kept in contact with the brass template that was used as a 
guide (Figure [2]-e). To make the narrow section of the coupon, a template 
with two grooves was used. After the first side was machined, the sample 
was rotated about one screw and set along the other groove to machine the 
other side (Figures [2j-f and g). During the entire procedure, the orientation 
of the wafer relative to the original bone was tracked and the bone was kept 
hydrated (0.9 % saline solution). The coupons were then stored in a tube 
filled with saline solution and kept in a fridge until the test day (between 
1 and 4 days after preparation) to ensure proper conservation, and imaged 
with the micro-tomograph. The thickness and width at the center of each 
coupon were measured prior to each test. 

2.2. Test fixture and procedure 

A hydraulic tensile machine was used (Model 8874, Instron Inc, Nor- 
wood, MA, USA). A clamping system made of aluminum was designed to 
avoid misalignment between the top and the bottom ends of the coupon 
(figure [4| . Because of their fragility, it was critical to be able to install the 
coupons on the test machine without applying any load that could gener- 
ate bending, shear or too much tension. Conventional fixed wedge clamps 
would not ensure the required control of the load during the installation of 
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the coupon. Therefore small low-mass clamps (10.9 grams each) were de- 
signed to be attached to the coupon (figure [5]) prior to mounting the sample 
to the test machine. A torque of 1 Nm was applied to the screws used to 
clamp the samples. Preliminary tests showed that this torque value was 
adequate to prevent slippage of the sample, while avoiding bone crushing. 




Figure 4: Clamping system for the coupons. 



Pins that go through the rod-end ball joints were used to connect the 
clamps to the slitted blocks affixed to the base of the tensile machine (bottom 
clamp) and to the end of the piston (top clamp, Figure [4]). 

The top clamp was first installed and the coupon equipped with the 
clamps was let to hang. The position of the sample was adjusted by moving 
the piston so that the bottom pin could go through the bottom ball joint 
without applying any load. The piston was then slowly moved up until 
the bottom pin touches the clevace, resulting in a preload of about 2 N. 
The bone coupons were stored in saline solution until the paint pattern was 



applied, a few minutes prior to testing (see section 2.4) 



Screw to apply the 
clamping force 




Figure 5: Close-up view of the clamp. 



2.3. Test matrix 

The sixth and seventh right ribs were harvested from the three sub- 
jects included in Lessley et al. (2010a) (Table [l]). Ten bone coupons were 



7 



machined, and tested up to fracture with a constant velocity displacement 
under quasi static (0.01 mm/s and 0.02 mm/s) or dynamic (24 mm/s) rates 
(Table [2]) . The velocity of the applied displacement was determined to gen- 
erate a strain rate similar to that reported on in Lessley et al. ( |2010a[ ) for 
the dynamic tests. 



Subject 


Age at time of death 


Cause of death 


Body mass (kg) 


Stature (cm) 


468 


67 


Stroke 


64 


166 


473 


54 


Brain aneurysm 


73 


182 


480 


71 


Laryngeal cancer 


70 


182 



Table 1: Subjects characteristics. 



Subject 


Rib level 


Location 


Aspect 


Velocity (mm/s) 


Strain rate (%/s) 


468 


6 


A 


Lateral 


0.02 


0.089 




7 


B 


Lateral 


0.01 


0.055 




7 


C 


Lateral 


24 


170 


473 


6 


B 


Medial 


0.01 


0.046 




7 


B 


Lateral 


0.01 


0.043 




7 


B 


Medial 


24 


170 




7 


C 


Lateral 


24 


170 


480 


6 


A 


Lateral 


0.01 


0.069 




7 


B 


Lateral 


24 


170 




7 


C 


Lateral 


0.01 


0.055 



Table 2: Test matrix. The actual strain rates are provided for the quasi-static tests, 
whereas only the target strain rate is provided for the dynamic tests as the strain time 



history data could not be determined (see section 2.4| 



2.4- Data acquisition and processing 

The tensile load was measured by a three axis loadcell located under the 
bottom clamp (model 6085, Denton Inc, Plymouth, MI, USA), connected to 
a standard data acquisition system (DEWE-2010, Dewetron GmbH, Graz, 
Austria). The tensile and shear loads were measured. The sampling rate 
was 500 Hz for the quasi-static tests, and 100 kHz for the dynamic tests, 
and the signals were filtered using a low-pass 2 nd Butterworth filter with 
a cut-off frequency of 10 Hz for the quasi-static tests and 600 Hz for the 
dynamic tests. These filters were chosen to remove the noise of the signals 
while keeping their shapes and introducing negligeable time-shift. 
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For the quasi-static tests, the strain was measured by performing image 
analysis of pictures taken during the tests with a 12-megapixel single- lens 
reflex camera equipped with a 100- mm macro lens and triggered by a pro- 
grammable controller (The Time Machine, Mumford Micro Systems, Santa 
Barbara, CA, USA) that allows for tripping the camera shutter at a constant 
frequency of 3 frames per second. The controller signal was also sampled by 
the data acquisition system to synchronize strain and force measurements. 
The outermost surface of the coupon (with respect to its orientation rela- 
tive to the rib before machining) was painted with a black and white pattern 
(figure [6]) to a llow for strain measurement based on contrast analysis (Frank 



and Spolenak , 2008 ) . Strain could not be measured for dynamic tests due to 



the unavailability of a continuous (non-flickering) light source. The tensile 
stress was defined as the ratio of the tensile force by the cross-sectional area 
in the middle of the coupon. The effective Young's modulus (referred to as 
Young's modulus) was defined as the slope of the strain-stress curve between 
and 0.5 %. For the quasi-static tests, stress and strain could be measured. 
For the dynamic tests, only failure stress could be determined. In addition, 
fracture location was documented by measuring the distance between the 
fracture line and the anterior end of the coupon. 



Figure 6: Coupon with high contrast pattern prior to test. 



2.5. Analysis of the images from micro- computed tomography 

The images obtained with the micro-tomograph (DICOM format) were 
analyzed to determine the actual cross-sectional area that accounts for the 
presence of voids in the cortical microstructure and for the cross-sections 
not being perfectly rectangular. A Matlab script (2010a, The Mathworks, 
Natick, MA, USA) was written to perform the segmentation of the DICOM 
images: a threshold of 3000 Housfield Unit (HU) was used, and all the pixels 
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with a HU value greater than the threshold were considered as bone. The 
position along the coupon was normalized, with and 1 being respectively 
the anterior and posterior aspects of the coupon (relative to the rib). 

3. Results 

3. 1 . Material properties 

Strain-stress curves were plotted for all the quasi-static tests (figure [7]). 
The material properties were determined by assuming a linear elastic model 
between and 0.5 % of tensile strain (table [3]). Young's modulus ranged 
from 11.4 to 18.5 GPa, failure stress from 83.4 to 143.9 MPa, and failure 
strain from 0.71 to 1.49 % in quasi-static, and failure stress ranged from 
94.7 to 155.9 MPa in dynamic. 




Strain (e) [%] 

Figure 7: Strain-stress curves for the quasi-static tests. *: sample fractured in the clamp. 

3.2. Fracture location 

Fracture occurred outside the grip area in nine of the ten tests. Fracture 
occurred first within the clamping area for one sample (473, 7, B) because 
the sample slid of the clamp, causing the screw to bear the tension load. The 
posterior end (right, figure [3]) for this sample was too thin (about 0.43 mm) 
to provide sufficient clamping load, and the sample slid relative to the clamp 
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Coupon 



E (GPa) a failure (MPa) e f ai i ure (%) Fracture location 



Quasi-static 



468, 6, A 


18.5 


129.9 


0.71 


468, 7, B 


14 


143.9 


1.19 


473, 6, B 


12.2 


83.4 


0.9 


473, 7, B 


12.8 


(*) 


(*) 


480, 6, A 


11.4 


106 


1.49 


480, 7, C 


12.2 


97.3 


1.02 


Mean values 


13.5 ± 2.6 


112.1 ± 24.5 


1.06 ± 0.29 




± standard deviation 



Dynamic 



468, 7, C 

473, 7, B 

473, 7, C 

480, 7, B 
Mean values 
± standard deviation 



140.8 

155.9 

107 

94.7 
124.6 ± 28.5 



Table 3: Rib cortical bone material properties for the quasi-static and dynamic tests, and 
images of the fractured coupons (posterior surface). The anterior end is the left side. (*) 
sample fractured in the clamp, (f) the fracture by the right hole occurred after the test 
during handling. 



after a certain tensile load was applied. The analysis of the images collected 
during the test for strain measurement confirms that sliding occurred after 
0.5 % deformation; therefore only the Young's modulus could be determined 
for this test. The variation of the cross-section along the length of the 
coupons shows that the gage area cross-sectional areas was close to the 
target (1.25 mm) for all the coupons (figures [8] and [9]) , whereas the cross- 
sectional area of the ends varied greatly from sample to sample. The images 
obtained from micro-computed tomography are provided for each coupon 
(Appendix A). 
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Figure 8: Variation of the cross-section (quasi-static tests). The fracture location is indi- 
cated by the dot. Zero is the anterior end of the coupon, and one is its posterior end. 



4. Discussion 

4-1. Material properties 

The rib bone mechanical properties were calculated based on the as- 
sumption that the rib cortical bone was homogeneous and isotropic. The 
average value and standard deviation for the Young's moduli determined for 
the quasi-static tests (13.5 ± 2.6 GPa) were similar to the 13.9 ± 3.7 GPa 



reported in Kemper et al. (2005). The failure stresses were also in the same 



range for both the quasi-static (121 ± 24.5 MPa) and dynamic tests (124.6 



± 28.5 MPa) in the current study compared to Kemper et al. (124.3 ± 35.4 
MPa). Although the target strain rates were different in these two studies 
(0.07 %/s for the quasi-static and 170 %/s for the dynamic tests in the cur- 



rent study, 50 %/s in Kemper et al. ), the measured strain rates are actually 
wide spread (0.043 %/s to 0.089 %/s for the current study, and from 9 %/s 
to 90 %/s in Kemper et al. ). Based on Hansen et al. (2008) analysis of sev- 



eral studies reporting bone material data, the Young's modulus is little rate 
sensitive for strain rates below 200 %/s. Therefore the comparison of the 
results reported in the current study with that reported in Kemper et al. is 



justified. A major difference between these two studies is the non linearity, 



as Kemper et all reported on a substantial yield behavior, and much greater 
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Figure 9: Variation of the cross-section (dynamic tests). The fracture location is indicated 
by the dot. Zero is the anterior end of the coupon, and 1 is its posterior end. 



failure strains (2.68 ± 1.4 %) than what was determined here (1.05 ± 0.29 
%). In addition, the failure strains reported on in Kemper et al.| range from 
0.52 to 6.56%/s. This range is high compared to what is commonly used for 
bone. It is not clear whether the discrepancies are artifacts because of the 
experimental methods or due to age differences, as the subjects in the cur- 
rent study are a little bit older than in Kemper et al. (2005). The coupons 
machined in Kemper et al. had a wide range of cross-sectional areas (from 
0.38 to 2.08 mm 2 ), whereas the cross-sectional areas of coupons used for the 
present tests were more uniform (1.25 mm 2 ). The effect of the size of the 
coupons on the measured mechanical properties is not well understood, in 
particular how the loading varies with the varying cross-sectional area. In 
the analysis performed in the two studies, tensile is assumed to be the main 
type of loading. Thanks to the image analysis performed to measure strain 
in the current study, it was possible to ensure that there was no slipping in 
the clamps (except for the one test identified in the results section) , whereas 



Kemper et al.l reported that most of the coupons broke in the clamp area 



which suggests that these fractures were artifactual of the experimental set- 
up. Too few samples were tested in this study to determine whether bone 
properties vary as a function of the location along the rib. 
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Study 



Bone type 



Young's modulus (GPa) Failure stress (MPa) Failure strain (%) 



Yamada and Evans| 
Yama da and Bvans| 
Hansen et all J20081 



1970[ ) 
M7| 



Kemper et al.| 1|?005 ) 

Current study - Quasi-static 
Current study - Dynamic 



Dry Human long bones 18.3 
Wet Human long bones 21.1 

Human femora 16.1 ± 2.1 

Human ribs 13.9 ± 3.7 

Human ribs 13.5 ± 2.6 
Human ribs 



140 

172 
119.8 ± 20.7 
124.3 ± 35.4 
112.1 ± 24.5 
124.6 ± 28.5 



1.49 
1.3 
2.5 ± 0.8 
2.68 ± 1.4 
1.06 ± 0.29 



Table 4: Summaries of the cortical bone material properties reported by various authors. 



The failure stress was found to be on average higher for the dynamic 
tests than for the quasi-static tests (however fewer samples were tested un- 
der dynamic loading). There is very little data in the published literature 
regarding rib cortical bone rate sensitivity under tensile loading. The data 
reported on in McElhaney ( 1966 ) are often mistakenly used to describe the 
rate sensitivity of the bone material properties in tension, although the re- 
sults established by McE lhaney are based on compression tests. Other stud- 
ies such as Crowninshield and Pope (1974) used animal tissue and the bone 
samples were machined in the lower extremity; consequently the coupons 
are bigger. In addition, the mechanical function of a a bovine leg is different 
from that of a human rib, and therefore the microstructure is likely to be 
different. It is indeed true that the coupons machined in the human cortical 
femur and tibia bone are denser than that machined in the rib bone (on- 
going analysis of unpublished test data). Therefore, the effective modulus 
is higher for the lower extremity coupons than for the rib ones (table [4]). 

The coupons machined in the human cortical bone are heterogeneous 
(Appendix A), and the actual cross-sections (determined by taking into 
account the presence of the voids) were smaller than the target of 1.25 mm 2 
(figures [8] and [9J . When compared to the material parameters determined 
for the human femur (table [4]), rib bone exhibits a lower Young's modulus, 
but the failure stress and failure strain are in the same range. Although 
the strains were not measured using the same methods for all these studies 
(extensometer or non contact optical measurements), it makes sense that 
the effective modulus appears smaller for the rib bones. The voids in the rib 
bone microstructure are not accounted for in the measurement of the cross- 
sectional area as the coupon cross-section is assumed to be rectangular and 
its area is approximated by measuring its thickness and width. With the 
femur cortical bone being denser than the rib cortical bone, coupons made 
of femur bone or rib bone with the same external dimensions would appear 
to have the same cross-sectional area, whereas the actual cross-section of 
the rib bone coupon will be smaller because of the voids. 
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The analysis performed in the present study relies on the assumption 
that bone is homogeneous and isotropic, and therefore that the stress is 
greater in the gage area. However, the fracture location was not systemati- 
cally in the gage area (table [3]), which suggest that the continuum mechanics 
approach has its limitations. This has several implications in terms of frac- 
ture prediction, as the microstructure - in addition to the geometry - alters 
the stress field in the cortical shell and therefore the fracture threshold and 



location. Li et al. (2010b) performed an extensive computational analysis 
to determine which features of the rib geometry and structure needed to 
be included in a finite element model in addition to the adequate material 
properties to predict fracture under antero-posterior loading. Load and dis- 
placement at fracture, could be successfully predicted, whereas the fracture 
location predicted by the model did not match the experimental results. 
Both the cortical and trabecular bones were assumed elastic homogeneous 
and isotropic, and this may be why the fracture location was not correct, 
even with the cortical thickness distribution mapped from the actual ribs 
used in the experiments. 

4-2. Fracture prediction 

None of the published studies on bone coupons reported the fracture 
locations. However this would provide a valuable piece of information re- 
garding the modeling approach that should be used to adequately predict 
bone fracture. As reported in the current study, fracture did not always took 
place in the gage area. Preliminary tests were performed on plastic homoge- 
neous and isotropic plastic coupons prepared following the same procedure, 
and they all fractured in the gage area because of tension (the fracture line 
was perpendicular to the long axis of the coupons). Besides, the shear forces 
recorded with the three axis load cell were small compared to the load mea- 
sured in the tensile direction (less than 5 %). This confirmed that the test 
apparatus itself did not generate artifactual fractures, and therefore the frac- 
ture characteristics are the results of the coupons material and structural 
properties. Therefore, although bone material might be adequately assumed 
to be homogeneous and elastic linear for sub- fracture behavior, the fracture 
properties seem to be dependent on the bone microstructure which is not 
properly accounted for with an elastic linear model. Researchers have had 



limited success to predict fracture location in finite element models ( Li et al 



2010b). Current finite element models of the ribs include the rib as structure 
composed of two distinct materials (trabecular and cortical bone). The re- 
sults of the current study suggest that the delineation between cortical and 
trabecular bone has become less clear with improving imaging capabilities 
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(figure [I]): the cortical shell in the rib has pores, and the transition from the 
cortical to the trabecular bone can be somewhat blurred. Rather than using 
a two-part model for bone, a gradient approach may prove more accurate to 
predict the actual deformation modes and fracture locations under dynamic 
loading. Rib bone material properties need to be refined to include a more 
precise model for fracture prediction, such as stress concentration caused by 
the presence of pores or by the connection between the trabeculae and the 
cortical shell. Micromechanics approaches such as the cohesive zone models 
( Subit et al. 2009 ) would be worth evaluating as they can predict the onset 
of a crack based on the local microstructure. 



4-3. Experimental set-up 

The test apparatus designed for this study is promising: by design no 
load is applied when the coupon is affixed to the machine, and only tension 
is applied during loading. This led to a very high success rate during testing 
(no coupons were broken during handling or connection to the machine). 
However, it is not clear how the "loose" boundary conditions (compared to 
the traditional wedge clamps) affect the fracture outcome. When a crack 
initiates, the load distribution is the coupons changed suddenly and the 
clamps are likely to reorient themselves, which could lead to premature 
fracture. The same phenomenon could occur with plasticity. None of the 
data collected in the current study allow for checking whether the clamps 
accelerate the fracture process. 

A non contact strain measurement system provided a non invasive method 
to estimate strain in the coupon. With the current method, only the av- 
erage strain was estimated, similar to what a strain gauge would measure. 
This prevents from determining the link between the microstructure and 
the strain field, and therefore the fracture mechanism. Finally, the strain 
measurement procedure is 2D only, and consequently any relative motion of 
one of the clamps towards or away from the camera was seen as a change in 
strain, and therefore the measured strain could be overestimated (and the 
Young's modulus underestimated). 



5. Conclusion 



Predicting rib fracture based on the computational model of the thorax 
remains an elusive challenge: although the contribution of the geometry to 
the fracture mechanism has been demonstrated and included to some extent 



in finite models of the thorax (Song et al. 2009 Li et al. 2011), the con- 
tribution of the bone material properties has to be better described. The 
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intrinsic bone features that contribute to the onset of fracture are not well 
comprised by the commonly used material models (isotropic and homoge- 
neous). The paradigm that describes bone as two entities — trabecular or 
cortical — needs to be revised to include a finer description of the bone 
microstructure. The results presented in this paper supplemented with past 
research highlight that the bone fracture mechanism is not well accounted 
for with a linear elastic model, at least for the rib bone. The attempt made 
in the current study to capture the bone microstructure by measuring the 
effective cross-sectional area along the coupons length proved unsuccessful. 
A more in-depth analysis of the bone microstructure (such as the direction 
of the voids and pores in the bone cortical shell) and how it modifies the 
strain field compared to the simplified approach that considers bone as elas- 
tic linear is required. It would allow to determine how the microstructure 
generates areas of weakness and strength to establish whether local strains 
or stresses play a role in the fracture mechanism. 
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Figure A. 10: Images of the cross-section obtained from micro-computed tomography 
(coupons tested under quasi-static loading) . Only one every ten images is included. 
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Figure A.ll: Images of the cross-section obtained from micro-computed tomography 
(coupons tested under dynamic loading). Only one every ten images is included. 
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